Mott transitions in correlated electron systems with orbital degrees of freedom 
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Mott metal-insulator transitions in an A/-fold orbitally degenerate Hubbard model are studied by 
means of a generalization of the linearized dynamical mean-field theory. The method allows for an 
efficient and reliable determination of the critical interaction U c for any integer filling n and different 
M at zero temperature. For half-filling a linear dependence of U c on M is found. Inclusion of the 
(full) Hund's rule exchange J results in a strong reduction of U c - The transition turns out to change 
qualitatively from continuous for J = to discontinuous for any finite J. This result is confirmed 
using the exact-diagonalization method for M — 2. 

PACS: n.lOFd, 7f.27.+a, 7f.30.+h 
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I. INTRODUCTION 

The Mott metal-insulator transition (MIT) has been 
the subject of numerous experimental and theoretical in- 
vestigations over the past decades 1-3 . As a transition 
from a paramagnetic metal to a paramagnetic insulator 
that is driven by the competition between the conduc- 
tion electrons' kinetic energy and their mutual Coulomb 
repulsion, the Mott MIT represents a prime example for 
a quantum phase transition. Realizations of the Mott 
transition can be found in different 3d transition-metal 
oxides 1 ' 2 as well as in alkali-doped Fullerides 4 ' 5 , for exam- 
ple. Its description is a challenging task for many-body 
theory. 

The minimal model to describe the Mott transition is 
the Hubbard model 6 which in its simplest form contains 
the local Coulomb interaction between the electrons in 
a conduction band formed by a single orbital only. The 
investigation of such a single-band model might be jus- 
tified in some cases, in particular as an effective model 
for low-energy properties. A consistent description of the 
experimental observations for the above-mentioned sys- 
tems, however, certainly requires the use of more realistic 
models which include orbital degrees of freedom, possible 
crystal field splittings, etc. 

Theoretical advances in the past decade - mainly due 
to the development of the dynamical mean-field the- 
ory (DMFT) 7,8 - have led to an increased understand- 
ing of multi-orbital Hubbard models. Recent investi- 
gations have either concentrated on fundamental ques- 
tions such as the nature of the Mott MIT as a func- 
tion of orbital degeneracy 9-16 or, within the so-called 
LDA+DMFT approach 17 , on a realistic description of 
the transition. 

Nevertheless, there are still no reliable DMFT results 
for the zero-temperature (T = 0) Mott transition in 
general multi-band models: Numerical methods to solve 
the DMFT equations are either restricted to fairly large 
temperatures (quantum Monte-Carlo approach, QMC) 
or have not yet been extended to the multi-band case 



(numerical renormalization-group method, NRG). Addi- 
tionally, due to the sign problem within the QMC ap- 
proach only a simplified exchange part of the Coulomb 
interaction can be taken into account 13,18 . The exact- 
diagonalization method (ED) is restricted to rather small 
system sizes, and data are available for M < 2 or- 
bitals only 16 . The M dependence of the critical inter- 
action strength U c for different integer fillings n is par- 
ticularly interesting. Detailed results for U c as a func- 
tion of M and n are available from the Gutzwiller ap- 
proximation and from a slave-boson approach 9,11 . Nu- 
merical DMFT-QMC calculations have been performed 
for M < 3 12,13 . Remarkably, in the limit of large or- 
bital degeneracy M — > oo an analytical treatment of 
the DMFT becomes possible for the MIT 15 . A scaling 
U c = £4,2 oc M for the actual transition is found while 
U c ,i v~M is obtained for the critical interaction where 
the insulating solution breaks down 15 . This is consistent 
with the linear dependence for large M found in Refs. 9, 
11 and with the square-root dependence reported in Ref. 
5. However, there is still a need for a DMFT method 
which works at T = and which allows for an efficient 
and reliable determination of U c for arbitrary n and M 
and for a model including the full exchange interaction. 

Here we present and apply a proper multi-band gener- 
alization of the so-called linearized DMFT (L-DMFT) 19 
which is a simple but rather successful technique that 
originally was developed for the critical regime of the 
Mott MIT in the single-band model. Within the L- 
DMFT the lattice problem is mapped onto an Anderson 
impurity model with a single bath site only by consid- 
ering a simplified self-consistency condition just at the 
critical point. The approach can be considered to be the 
simplest non-trivial variant of the more systematic pro- 
jective self-consistent method (PSCM) by Moeller et al 20 . 
It allows for extremely fast numerical calculations or even 
analytical results to estimate the critical parameters. Ex- 
tensions have been considered for the (Mott-Hubbard or 
charge-transfer) MITs in a d-p model 21 , for a quantum 
critical point in the periodic Anderson model 22 , and for 
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the M = 2 degenerate Hubbard model at half-filling and 
J = 16 . Furthermore, a general two-site DMFT 23 has 
been developed which reduces to the L-DMFT at the 
MIT but is not restricted to the critical point. An appli- 
cation to multi-band systems is likewise possible but has 
not yet been considered. 

As concerns effective single-band models, the linearized 
DMFT has been tested extensively by comparing with 
numerical results from the full DMFT. For the standard 
Hubbard model 19 ' 24 but also for different thin-film and 
semi-infinite surface geometries 25 as well as for the d-p 
model 26 a remarkable agreement has been found: De- 
tailed trends of U c as a function of electronic model pa- 
rameters and as a function of parameters characterizing 
the lattice geometry are predicted reliably. For quanti- 
tative estimates errors of the order of a few percent have 
to be tolerated. 

As it is demonstrated here, a straightforward extension 
of the method to multi-band systems is possible. We ex- 
pect the L-DMFT to be a fast but reliable tool in this 
case as well since the effective impurity model includes 
the complete atomic part and the self-consistently de- 
termined bath with the complete degeneracy (although 
restricted to a single site only). Furthermore, the relia- 
bility of the L-DMFT for the multi-band case is checked 
by comparing with different fully numerical DMFT tech- 
niques whenever results are available. All DMFT calcu- 
lations discussed in this paper were performed on a Bcthc 
lattice with t 2 = 1, corresponding to bandwidth of 4. 



II. MODEL AND LINEARIZED DMFT 

We consider a multi-band Hubbard-type model H = 
H + Hi including an intra-orbital hopping 

H = ^j^L/j™ (!) 

ijOLO 

and the direct and the exchange part of the on-site 
Coulomb interaction: 

H\ = — U aa ' c\ a(y c\ aliy ,Ci a r a iCi aa 

+ ^ '^ aa ' C iaa C ia'a' C ia<T' c ia'(T ; (2) 

iaa' ' aa' 

using standard notations, i is a site index, and a =T,I 
refers to the spin direction. The different orbitals la- 
beled by a = l,...,M are considered to be equivalent, 
M is the orbital degeneracy. Exploiting atomic symme- 
tries, we have U aa > = (U + 2J)S aa > + U(l — S aa >) and 
Jaa' — JO- —Saa') as usual (see e.g. Ref. 27). Opposed to 
calculations using the QMC method (see, e.g. Refs. 13, 
18), Hi includes the full exchange part and thus preserves 
rotational invariancc. 

The DMFT assumes the self-energy to be local. Fur- 
thermore, the self-energy as well as the Green function 



are diagonal with respect to the orbital index a for the 
model considered here. In case of equivalent orbitals and 
in the absence of any spontaneous symmetry breaking, 
we thus have £ij a/ 3(a->) = 5ij5 a pY<(uj) and likewise for the 
local Green function: Gu a f3(uj) = S a pG(uj). Within the 
framework of the DMFT, the model H is mapped onto 
an impurity model H' = H' + H[ where H[ is the local 
part of the interaction (2) at a distinguished site io and 

n s 

H = X! to C ioaa C toaa + ^ £fc a kaa a kaa 
cut aa.k=2 

"a 

+ E V k (c\ aaa a kaa +\i.c.) (3) 

acr,fe=2 

with to — tu. The self-energy £(w) of H' is identified 
with the self-energy of the lattice model H and yields, 
via the lattice Dyson equation, the on-site Green func- 
tion G(uj). The latter determines the one-particle pa- 
rameters of H' or, equivalently, the hybridization func- 
tion A(w) = J^k V k l{u + M - £fe) via the DMFT self- 
consistency condition A(w) = uu + n — to — ^uj) — 1/G(u>). 

Within the full DMFT an infinite number of bath de- 
grees of freedom n s oo are necessary to fulfill the self- 
consistency equation. This implies the need for further 
approximations to treat H' . Contrary, a simple two-site 
model is considered here. For n s — 2 the model (3) can 
be solved exactly; however, a simplified self-consistency 
condition must be tolerated. 

Our approach rests on the assumption that similar as 
for the single-band model, the MIT is characterized by 
a vanishing weight of a low-energy quasi-particle reso- 
nance. Close to the MIT it is then plausible to neglect 
any internal structure of the quasi-particle peak and to 
assume that there is no influence of high-energy features 
on the low-energy part of the excitation spectrum. The 
quasi-particle peak is then approximated by a single pole 
at the Fermi energy (e c = e k =2 = m) which reproduces 
itself in the DMFT self-consistency cycle. One can pro- 
ceed as for the single-band model 19 and is finally left with 
the following simple algebraic self-consistency equation 
which determines the critical parameters (V = Vk=2)'- 

V 2 = zMf ] . (4) 

Here is the second moment of the non-interacting 

density of states = / dxx 2 p^{x). Formally, Eq. 

(4) is the same as in the single-band model 19 . However, 
the quasi-particle weight z has to determined from the 
degenerate model (3). We calculate z = 1/(1 — S'(0)) 
from the spectral representation of G(uj) by summing 
the residua oc V 2 in the limit V — > 0. The residua are 
obtained by standard (Brillouin-Wigner) degenerate per- 
turbation theory up to the order V 2 . This yields an ex- 
pression z = z{V, e c , t , U, J) = V 2 F(e c -t , U, J)/M 2 (0) + 
0(V 4 ). Using e c = [i in the limit V — > and setting 
t = tu = to fix the energy zero, we obtain from Eq. 
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(4) F(/i, U, J) = 1 as a condition for the MIT. This con- 
dition is different for different orbital degeneracy M and 
integer filling n — 1, 2 M — 1. 



III. RESULTS AND DISCUSSION 

For an integer filling n and for any U larger than a 
critical interaction U c = U c (n, M) there are two critical 
values for the chemical potential /J.±(U) which can be 
obtained by solving F([i, U) — 1 for \x (the case J = is 
considered first). For U > U c {n,M) and H-(U) < fi < 
H+{U) the system is a Mott insulator with integer filling 
n. On the other hand, for ^ < H-(U) or \i > ^i + (U), 
a non-integer filling is realized, and the system becomes 
metallic. The chemical potential is a continuous function 
of the filling n' in the vicinity of n for any U < U c (n 7 M) 
but n(n') is discontinuous at n' = n with a jump from 
H-{U) to n+{U) for U > U c (n,M). U c is determined by 
minimization of the function U (/x) which is obtained by 
solving F(fi, U) = 1 for U. 

For half-filling n = M and arbitrary M we succeeded 
with a completely analytical calculation for F(fi,U): 



F(»,U)/M. 



(0) = 1 ( M+l M_ 

2\MU-n n-{M-l)U 



M 



+ 



M + l 



2 \MU - ii n-{M-l)U 



(5) 



For U < U c and half-filling n = M in the symmetric 
model, the chemical potential is fixed to /i = (M — \)U 
due to particle-hole symmetry. Then we find F = 
M 2 (0) (4M + 2) 2 /U 2 which results in the critical interac- 
tion 



U c {n = M,M) = (4M + 2)^M. 



(0) 



(6) 



By solving U) = 1 with Eq. (5) for [/ > U c we 

obtain the jump in the chemical potential 



error due to finite-size effects in the ED is almost negli- 
gible. One can state that the agreement of the L-DMFT 
results with the ED and PSCM data is equally good for 
both M = 1 and M = 2. 

The generalized L-DMFT predicts U c (n = M, M) to 
depend linearly on M, see Eq. (6). This agrees with 
previous results from the Gutzwiller method 9 , slave bo- 
son calculations 11 , the non-crossing approximation for 
low temperature 28 , and the projective technique for large 
M 15 . This finding can easily be understood by looking 
at the weight factors in the Lehmann representation of 
the spectral density of the two-site Anderson model. Ex- 
panding into the different configurations, one finds that 
the ground state is composed of the order of M singlet 
states in the limit V — > and that there are of the order 
of M excited states coupling to the ground state. Hence, 
there are of the order of M 2 processes to be considered in 
the t^-perturbation theory which are contributing to the 
spectral weight near w = 0. This implies z oc M 2 V 2 /U 2 
and eventually results in U c cx M. 

In Fig. 2 the dependence of U c on both n and M is 
shown for all integer fillings and M < 5. We find U c to be 
reduced away from half-filling, consistent with previous 
results 5 ' 9,12-14 . The fit shows an almost perfect quadratic 
dependence of U c on the filling. 

Let us now discuss the influence of the Hund's rule 
coupling J on the Mott transition at T = 0. Contrary 
to the QMC method, the full exchange part of the on- 
site Coulomb interaction can easily be included in the 
L-DMFT calculation. A quantitative influence, i.e. a re- 
duction of U c due to J has already been observed in Refs. 
5, 13. This behavior is also evident from our L-DMFT 
results, as shown below. 

The situation is, however, more complicated: It turns 
out that a finite J leads to a qualitative change of the 
nature of the T = Mott transition. This is illustrated 
in Fig. 3 which shows the [/-dependence of the quasi- 
particle weight z for various values of J and n — M = 2. 
Note that z is plotted not only close to the transition 
but also for U far below U c which means that the self- 



A M ^ + - M _ = [/yB-yi? 2 -l + p (7) 

where B = 1 + 1/(2(2M + \) 2 {U/U C ) 2 ). The result for 
A/j,(U) is shown in Fig. 1 for M — 1,...,5. Here and 
hereafter we set — 1. For fixed M and U — > oo, the 
jump in the chemical potential increases linearly with U. 
For U — ► U c , A/i(U) shows a square-root behaviour. 

Qualitatively and even quantitatively, this agrees well 
with results from numerical methods to solve the full 
DMFT equations. In particular, we compare with re- 
sults from ED calculations which we have performed 
with n s — 6 for M = 2 (circles) and with data avail- 
able for T = from the projective self-consistent method 
(crosses) 14 . The non-degenerate case M = 1 has already 
been discussed in Ref. 24. For M = 2 and n s = 6 the 
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FIG. 1. Afi(U) for half-filling n = M and different M as 
obtained from Eq. (7). Results from the full DMFT using 
ED, NRG, and PSCM 14 ' 20 are shown for comparison. 
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FIG. 2. Critical interaction U c for different fillings 
n = 1,2, ...,2M — 1 and different M as a function of 
n/M. Symbols: L-DMFT results. Lines: simple fit func- 
tion U c (n,M) = U c (n = M,M) — c(l - (n/M)) 2 and 
c= -0.217 + 2. 671M. 




FIG. 3. [/-dependence of the quasi-particle weight 2 for 
n = M = 2 and different J. Dashed lines correspond to 
metastable solutions; vertical lines indicate the actual Mott 
transition at U = U c . 



consistency relation (4) is used for any U. This is rea- 
sonable as long as there is a clear separation between 
the low-energy (quasi-particle) and the high-energy scale 
(Hubbard bands). Although the approximation becomes 
questionable for higher z and particularly in the limit 
z — > 1, the trend of z(U) looks rather plausible in whole 
range. 

In the critical regime where z goes to zero continu- 
ously, we find z(U) to decrease with decreasing U which 
seems to be unphysical. In fact, this behavior belongs 
to a metastable state (dashed line) while for the true 
ground state (solid line) z(U) is always decreasing with 
increasing U. The two metallic states with finite z > 
but different slope dz / dU are coexisting in a finite range 
of interaction strengths, namely for U larger than U c _2 
(given by z = 0) and for U smaller than another critical 
value ?7 Cj 3 > U c ^ where they merge together 29 . 
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FIG. 4. Energy of the two metallic solutions relative to the 
energy of the insulator for J = 0.1. At U c ,2 (closed circle) 
the metastable metallic state continuously merges with the 
insulating state. At U c , 3 (cross) the two metallic states merge. 
The discontinuous phase transition takes place at U c (open 
circle) . 
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FIG. 5. Local moment at the impurity site (S 2 ) as a func- 
tion of U for different J. Dashed lines: metastable solutions; 
vertical lines: actual Mott transition for U = U c - 



The actual critical interaction strength U c where the 
transition to the insulator takes place must be deter- 
mined by comparing the respective energies with the en- 
ergy of the insulator. This is shown in Fig. 4 for J = 0.1. 
For any J > it is found that U c , 2 < U c < [/ Cj3 . At 
U c the energies of the metallic and insulating states are 
crossing as functions of U. Consequently, the transition 
at T = is discontinuous with a jump in physical quanti- 
ties like the double occupancy (n-\n{) = dE/dU. Equiv- 
alently, one can look at the local magnetic moment at 
the impurity site (S 2 ) (see Fig. 5). The local moment 
increases with U, typical for a strongly correlated nor- 
mal metal. At U c it jumps to (S 2 ) = 2 for any finite 
J. Similarly, the quasi-particle weigth z decreases with 
increasing U and drops to zero at U c from a finite value 
z > 0. In Fig. 3 this is indicated by solid vertical lines. 

The phase diagram in the (U, J)-plane is shown in Fig. 
6. The critical interaction for the actual transition is 
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FIG. 6. (U, J) phase diagram for n = M = 2. U c : criti- 
cal interaction for the MIT. U c ,3'- maximum U up to which 
a metallic solution is found. U c ,2'- critical U where z — > 
(U c ,2 = 14/3 - 3 J for J > and M< 0) = 1). 
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FIG. 7. DMFT-ED results (n s = 6) for the 17-dependence 
of the quasi-particle weight (left scale) and the local moment 
(right scale). J = and J = 0.1 for n = M — 2. 

plotted together with U Ct 2 and t/ Cj 3. U c is found to be 
very close to U c ,3 for all J- values. Note that already a 
moderate value of J ss 0.3 leads to a drastic reduction of 
U c from U c = 10 to U c w 5. 

Interestingly, similar results have been obtained be- 
forehand by a rather different method: Applying a gen- 
eralized Gutzwiller approximation to a degenerate two- 
band model, Biinemann et al. 10 found the Mott transi- 
tion to be discontinuous for any finite J > and to be 
continuous for J = only. While the Gutzwiller method 
(for J = 0) significantly overestimates the critical inter- 
action strength as compared to the L-DMFT, the reduc- 
tion of U c by a finite J is stronger. Common to both 
methods, however, is the qualitative change of the Mott 
transition for a finite J. 

We have additionally checked the nature of the Mott 
transition by a fully numerical evaluation of the DMFT 
equations using the ED method with n s = 6 sites for 
n = M = 2. The results for J = and J = 0.1 are 
shown in Fig. 7. The critical interaction U c is found to 
be somewhat smaller as compared to the L-DMFT result 



as shown in Fig. 6. More important, however, the ED 
calculations show the same qualitative trend: The slope 
of z(U) approaches — oo at a critical value U = U c ^. 
This clearly indicates a second metastable metallic solu- 
tion and a discontinuous transition and thus corroborates 
our results discussed above. Fig. 7 also shows the U de- 
pendence of the local magnetic moment where the effect 
is even more pronounced. Note that in the insulator (S 2 ) 
is close but not equal to the prediction of the linearized 
DMFT (S 2 ) = 2 (cf. Fig. 5). This is due to residual local 
fluctuations which are neglected by the L-DMFT for the 
insulating phase. 

The limit J — > which appears to be somewhat ex- 
ceptional can be analyzed easily. From Fig. 3 it is obvi- 
ous that the critical interaction U c ,2 at which z — > in 
the metastable phase shows a discontinuous jump from a 
value U Ct 2 w 5 < f7 Ci3 for J > to U c . 2 = 10 = C/ C)3 for 
J = 0. To find the physical reason for this discontinuity 
of £/ c .2, one can again look at the Lchmann representation 
of the spectral density for V — > 0. The main difference 
between the case J = and the case J > consists in 
the fact that for J > only those configurations with 
a triplet on the impurity site contribute to the (singlet) 
ground state in the limit V — > (see Fig. 5). As com- 
pared with J = this is about a factor of 2 configurations 
less. Thus, the discontinuous jump of Z7 C; 2 is caused by 
a strong suppression of orbital fluctuations due to a fi- 
nite J (a similar argument in a slightly different context 
has been put forward in Ref. 16). On the other hand, 
it should be stressed that the actual critical interaction 
U c remains continuous in the limit J — > as U c — > £4.3. 
Note that to uncover this subtle distinction it has been 
necessary to extend the theory beyond the z — > limit. 

IV. CONCLUSION 

To summarize, we have presented an extension of the 
linearized dynamical mean-field theory (L-DMFT) to 
multi-orbital Hubbard models including the full exchange 
term. This allows a reliable calculation of the critical 
interaction strength U c for the zero-temperature Mott 
transition analytically for all M at half-filling and nu- 
merically for all integer fillings n 7^ M and M < 5. We 
find a linear dependence of U c on the number of orbitals. 
Remarkably, a qualitative change of the nature of the 
T = Mott transition is observed for any finite J: the 
transition becomes discontinuous with a finite jump of 
the quasi-particle weight at U c close to ?7 Cj 3. 

For finite temperatures one might speculate that a 
finite J increases the tendency for a first-order Mott- 
transition. This could be of experimental relevance for 
the Mott-transition in transition-metal oxides. 

The L-DMFT for correlated electron systems with or- 
bital degrees of freedom is a handy and trustable ap- 
proach with a large potential for future applications. One 
might consider, for instance, the multi-band Hubbard 
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model investigated in the context of the LDA+DMFT 
approach for the transition-metal oxide V2O3 18 . The L- 
DMFT can be easily extended to answer questions about 
the importance of the eg-bands (neglected in Ref. 18) on 
the Mott transition, or the influence of the full exchange 
term. We have already observed a significant quanti- 
tative change of U c when excluding the spin-dependent 
part from the exchange term in (2). This should have 
consequences for a realistic description of the MIT in 
transition-metal oxides which are worth to be studied. 
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